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Abstract 

We study radiative corrections to 1? —> Xcii^i decays assuming the power counting 
~ for the charm-quark mass. Concentrating on the shape-function 

region, we use effective field-theory methods to calculate the hadronic tensor at 
NLO accuracy. From this we deduce a shape-function independent relation be¬ 
tween partially integrated B —> Xcii>£ and B —> Xuii'i spectra to leading power in 
llrrih, including first-order corrections in the strong coupling constant. This may 
provide an independent cross-check on the determination of the CKM element 

\Vub\. 



1 Introduction 


A central goal of the B physics program is to accurately determine the CKM parameter 
\Vub\- A complication is that experiments cannot measure the total rate for inclusive 
B — Xuii^i decays, because part of the available phase space is dominated by a much 
larger background from B Xcii^e decays. In fact, data for inclusive b ^ u transitions 
is available only in the shape-function region, where the hnal-state hadronic jet carries a 
large energy on the order of the 6-quark mass mb, but a relatively small invariant mass 
squared on the order of Agco^^b- 

The study of inclusive B decays in the shape-function region using soft-collinear effec¬ 
tive theory (SCET) has received much attention in recent years PQ1213IH El • Predictions 
for decay distributions are available in the form of factorization formulas which separate 
the physics from the three scales mb ^ V^A^cd^ S> Aqcd- At leading order in l/mt, 
the factorization formula takes the form 


H-J(^S. (1) 

The hard function H and the jet function J are perturbatively calculable functions 
depending on quantities at the hard scale mb and the hard-collinear (jet) scale -^/AgcD^b 
respectively. The shape function S' is a non-perturbative function defined in terms of a 
non-local HQET matrix element [H]- There are two basic strategies for reducing shape- 
function related hadronic uncertainties in the measurement of \ Vub\- The hrst is to extract 
the shape function in one process and use it as input for other processes, the second is 
to construct shape-function independent relations between different decay distributions. 
Common implementations of these strategies use B —^ in combination with B — > 

Xuii^e decay spectra [3 013 [13 [mini 113 . 

Assuming the power counting ~ AqcD’^b for the charm-quark mass, parts of 
the phase space for B —»• X^^vi decays lie in the shape-function region mi uni- Work 
performed in mi showed that the singly differential spectrum in a certain kinematic 
variable has much in common with the P+ spectrum in P —> Xuii^e decays. In particular, 
at tree level and excluding power corrections, this spectrum is directly proportional to 
the leading-order shape function. This raised the possibility of using data from inclusive 
decays into charm quarks to learn about the leading-order shape function. The analysis 
in mi concentrated on the classification of sub-leading effects in the Aqcd/ expansion 
at tree level, while the question of radiative corrections was left open. 

In this paper we calculate the perturbative corrections to P ^ Xcii^i decays in the 
shape-function region. We show that our one-loop result for the hadronic tensor can be 
written in the factorized form (0). Moreover, the hard function H and the shape function 
S are identical for inclusive b ^ u and b ^ c transitions; the charm-quark mass affects 
only the jet function J. This allows us to construct a simple, shape-function independent 
relation between the P —> Xcii^e and P —> decay spectra, which may provide an 

independent cross-check for the determination of |I4fe|- 

The paper is organized as follows. In Section |2 we discuss some aspects of SCET 
needed in our analysis. We use this to calculate the hadronic tensor at one loop in 
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Section El In Section E) we present results for the partially integrated spectrum needed 


in our phenomenological discussion and examine some issues related to the dehnition of 
the charm-quark mass. A relation between partially integrated 6 —>■ u and 6 —*• c spectra 
is derived and studied in Section El We conclude in Section El 


2 SCET for B XcB'i transitions 


In this section we review some aspects of SCET na im ca EO] needed to describe 
inclusive b ^ c transitions in the shape-function region. The effective theory facilitates 
the separation of scales and sets up a systematic expansion in the small parameter 



( 2 ) 


At the level of Feynman diagrams, this separation of scales is achieved by evaluating QCD 
integrals using the method of regions j2I], and the construction of SCET is closely related 
to this diagrammatic analysis. To apply this method one hrst identihes the momentum 
regions which give rise to leading-order on-shell singularities in loop diagrams. The 
integrand is expanded in A as appropriate for the particular region before performing 
the integral. Once all the regions are identihed, their sum is equal to the full theory 
integral, up to higher-order terms in A. 

Applying the method of regions to inclusive 6 —> u transitions in the shape-function 
region, where the jet momentum and the jet energy satisfy ~ m^A^ and E ~ m^, one 
hnds contributions from hard, hard-collinear, and soft regions. SCET is constructed in 
such a way that the hard-collinear and soft regions are contained in effective theory helds 
and operators, while the hard region is contained in Wilson coefficients multiplying these 
operators. For the 6 —> c transitions dealt with in this paper, we will always work in the 
kinematic region where the jet momentum and the jet energy satisfy r-u r-u rnl\^ 

and E ~ rrih. It is apparent that the set of regions is identical to that in the charmless 
case; one must replace p^ — ml in the hard-collinear propagators, but the A 

expansion, and thus the regions calculation, works the same. Therefore, the relevant 
version of SCET is very similar to that for charmless decays. The objects of interest are 
the SCET Lagrangian and currents, which we now discuss in turn. 

2.1 SCET Lagrangian and mass renormalization 

The leading-order SCET Lagrangian for a hard-collinear quark with mass me interacting 
with soft and hard-collinear gluons is (see for instance 1221 ESI uni) 



( 3 ) 


Here is a hard-collinear quark held, and the covariant derivatives are dehned as irt-D = 
in_d + gn-Ai,c + gn-Ag and zDhc = id + gA\^c- We have introduced two light-like vectors 
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n±, which satisfy n+n_ = 2. The Yang-Mills Lagrangian Lym for the soft (hard-collinear) 
sector is the same as in QCD, bnt restricted to soft (hard-collinear) fields. 

In massless SCET, the Lagrangian is not renormalized, in the sense that no new 
operators or non-trivial Wilson coefficients are indnced by radiative corrections. The 
reasoning for this was given in 1201, and involves showing that certain momentnm regions 
give rise to scaleless integrals. These argnments also apply to the SCET Lagrangian ©, 
becanse the A expansion is nnaffected by the presence of a qnark mass in the hard- 
collinear propagators, as we have emphasized above. 

On the other hand, mass renormalization plays a non-trivial role in onr analysis, and 
will be needed in the next section when we calcnlate the one-loop jet fnnction. We panse 
here to discnss mass renormalization in SCET. Later on we will stndy the differential 
spectrnm in the variable 

u = n_p — ml/n+p , 


where rric may be taken as the pole mass (in the massless case u rednces to the variable 

= U-p). In Section Fd.ll we will discnss alternative mass definitions which indnce a 
change in the jet fnnction of order as- 

Mass renormalization in SCET is closely related to the usnal QCD prescription, which 
follows from the observation that the self-energy diagram in SCET can be obtained from 
the A expansion of the corresponding QCD diagram. This has been pointed ont for the 
massless case in na, and for the massive case with Agcomfe in j21|. We have 

confirmed that it also holds for the case ~ Agcnm;,. In fnll QCD, the one-loop 
fermion propagator is 


G{p) 


I 

^-nic- E(p) ’ 


( 4 ) 


where the fermion self-energy reads 


S(p) =^Sy(pQ +mcE5(pQ. 


( 5 ) 


Analogonsly, the one-loop fermion propagator in SCET is 


Gdv) 


i 

u — T,^{u,n+p) 2 


( 6 ) 


For simplicity we consider a frame where p± = 0, snch that u = ri-p — rril/n+p. We 
obtain the SCET fermion self-energy Sg(M,n+p) by expanding the QCD propagator (jlj) 
to leading order in A and matching it with the SCET propagator ®, which gives the 
resnlt 

9 

777 

S5(n,n+p) = nSy(p2) + ^2 (Sy(/) + S5(/)) . (7) 

n+p 

Taking into acconnt mass renormalization, the renormalized fermion propagator in SCET 

is 


GdP) 


i 

u — Ej(m, n+p) 


_^ 

2 ’ 

n+p 


( 8 ) 
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where ^{ml) = 2mc5mc- The propagator has a pole for <(=> m = 0, from which 

we get 


n+p n+p 


2 (EK(ml) + Es(m2)) 




n+p Att \e \ J 3 
as the corresponding mass counterterm in the pole scheme. 


( 9 ) 


2.2 SCET transition current 

Unlike the Lagrangian, the SCET representation of the weak transition current involves 
non-trivial hard matching coefficients. The matching onto SCET takes the form UHlEn] 

^ /• _ 

- 75)6(0;) ds Ci{s,mb){^W){x + sn+)T^ h^X-) 

i=i 

3 

= ^ Ci{n+p,mb){^W){x)T>^ K{x-) ( 10 ) 

i=l 

where is the heavy-quark held dehned in HQET, fU is a hard-collinear Wilson line, 
and p is the momentum of the hard-collinear quark. The Dirac structures are chosen as 

71^ 

rr = { 7"(1 - 75 ), ^^"(1 + 75 ), ^(1 + 75 )}. ( 11 ) 

n_v 

One calculates the hard coefficients Ci by matching the one-loop corrections to the 
current from QCD to SCET. The QCD diagrams receive contributions from hard, hard- 
collinear and soft momentum regions. Since SCET is constructed to reproduce the 
results for the hard-collinear and soft regions, it is only the hard region of the QCD 
diagrams which contributes to the matching conditions. However, the Taylor-expanded 
integrand for the hard region does not depend on the hard-collinear scale so the 
matching conditions are the same as in the massless case. We can therefore read off the 
result for the coefficients Ci from The matching conditions also involve a current 
renormalization factor, which accounts for the divergent part of the hard diagrams. Its 
explicit form is [TH] 


Zj = 1 + 


CpOts 

dvr 


- A 

e ^ /i 2e 


( 12 ) 


We will need this renormalization factor in our calculation of the hadronic tensor in the 
next section. 
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3 Hadronic tensor at one loop 


In this section we calculate the one-loop corrections to the hadronic tensor for B —> 
decays in the shape-function region, always working to leading order in A. The hadronic 
tensor contains all the QCD effects in the semi-leptonic decay and is the starting point 
for deriving differential decay distributions. We define the hadronic tensor as 

= -\m{B{v)\T^^'^\B{v)), (13) 

71 

where we use the state normalization {B{v)\B{v)) = 1. The current correlator is 
given by 

(14) 

where = C 7 ^(l — 75)6 is the flavor-changing weak transition current discussed above. 
The one-loop result for the hadronic tensor can be written in the factorized form 






*,j=i 


Hij{n+p) / dujJ{u — uj,n^p)S{uj) 


(15) 


where p = m^v — q is the jet momentum in the parton model. The hard functions Hij, the 
jet function J, and the shape function S contain physics at the scales Agco^fe, and 
Aqcd; respectively. The limits of integration in the convolution integral are determined 
by the facts that the shape function has support for —A < a; < cx) and the jet function 
has support for m — a; > 0 . 

The procedure leading to m is familiar from charmless decay and involves a two- 
step matching procedure 1110 . In the first step, one integrates out hard fluctuations at 
the scale mb by matching the hadronic tensor calculated in QCD onto that calculated 
in SCET. The associated matching coefficients are the hard functions Hij. Since these 
coefficients take into account the hard region of the QCD diagrams, and this region is 
unaffected by the presence of a quark mass in the hard-collinear propagators, they are 
identical to those in the massless case. One finds Hij = CjCi, where the Ci are the 
hard Wilson coefficients defined in m- In the second step, one integrates out hard- 
collinear fluctuations at the scale Agco’^fe by matching the hadronic tensor calculated 
in SCET onto that calculated in HQET. The matching coefficient from this step is the 
jet function J. This function is obviously more complicated than in massless SCET, 
since it can depend on as well as We will calculate it in the following subsection. 
However, the final low-energy theory is still HQET, and the matrix element defining the 
shape function is the same as in charmless decays. For this reason, we can write our 
result in the form m- 


3.1 One-loop jet function 

The calculation of the one-loop jet function is conceptually identical to that for the 
massless case 00 , and we will closely follow the treatment in [2j. The jet function 
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Figure 1: The one-loop SCET graphs contributing to the current correlator. Mirror 
graphs are not shown. Graph (d) shows the insertion of a counterterm from mass renor¬ 
malization. 


is the matching coefficient between the hadronic tensor calculated in SCET and that 
calculated in HQET. The relevant SCET diagrams are shown in Fig^ We calculate them 
in the parton model, using on-shell heavy quark states carrying a residual momentum k 
satisfying u/c = 0. We work with dimensional regularization in d = 4 — 2e dimensions, 
using the Feynman gauge. The result for the graphs involving hard-collinear gluon 
exchange, including the counterterm from mass renormalization in the pole scheme o, 
can be written as 


D 


( 1 ) 

he 


y(l) 

'-'he 




(16) 


where 


(1) _ CpO-s * j 4 3 4 


q-'y'-) _ 

'-^hc ~ 


, „ H-In 

47r u' 1 e e 

^2 


TT 


+ 7---31U 


n+p u 

/i2 


-n+pu 


+ 2 In^ 


-n+pu 


2 TT^ 


nir 


4Li2 1 + ^ 


w 


rUp nip {nip + 2 u') 


m, 


+ u' {nip + u'Y 


In 


u 


nir 


(17) 


Here u' = u + n_k, as = «s(p), and nip = ml/n+p. We have checked that our result (fT7j) 
agrees with the corresponding result in pg when expanded in and translated to 
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the MS scheme. The graphs involving soft gluon exchange give 




K Tj 




2 ^ I "'V 


where 


(1) ^ Cpas i 
47r u' 


2 2 4, 

-^ H-1— In 

e e 




(18) 


(19) 


The sum of the 1/e poles in is removed by current renormalization in SCET, 

which is implemented by applying a factor of Zj (see IT^ to the bare current correlator. 
This renormalization factor is related to the divergent part of the hard region of the QCD 
diagrams, which was integrated out in the hrst step of matching. That it cancels the 1/e 
poles from the SCET diagrams, which are due to both hard-collinear and soft regions, 
shows that we have indeed constructed the appropriate version of SCET. Moreover, the 
pole structure for each individual region is the same as in the massless case. It follows 
that the hard and shape functions obey the same renormalization group evolution as in 
the massless case, a fact which we will use when discussing decay distributions in the 
next section. 

We can interpret the imaginary part of the hnite pieces of the SCET diagrams as 
one-loop corrections to the factorized expression They take the form ® S^^rt + 
<8) where the superscript {n) denotes the n-loop correction to each function, 

and the ® stands for a convolution. The tree-level functions are = 5{u — oj) and 
S^art = 5{ijj+n-k). As in the massless case, the one-loop correction to the shape function 
in the parton model is related to To show this, we take its imaginary part, which 

can be expressed in terms of star distributions, dehned as 


oM 


duF{u) 


'<0 



, F{u) - F(0) 

du - 

u 


'0 


-I- F(0) In 



( 20 ) 


7^/ \ /ln(M/m) 
duF{u] ' ^ > 


'<0 V * Jo 

It is not difficult to derive the following formulas 


rT T(“)-^(o) in"+Ahi„Tj£y 

0 u m2 \m J 


( 21 ) 


1 . 


-Im 

71 

Im 

TT 


In 


u \ 1 
mJ u 

In^ ( - “ i 1 


m/ u 


= 2 


ln(M/m) 


u 


- y<^(«), 


( 22 ) 
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where to take the correct branch of the logarithms we reinstored u = u + ie. We then 
obtain 


j'“> ® = SS<(«') = -Im 


( 1 ) 


CpCts 

An 


5{u') + 4 


TT 


U' 




( 1 ) 


s, finite 


[m] 


+ ' 


Hu' //i) \ 


M' 


w 




(23) 


which is identical to the one-loop result calculated in HQET. The jet function is related 
to the imaginary part of the hnite piece of j2]. In particular, we have 


^ 71 


ijp 


( 1 ) 


hc^finite 


(24) 


Using (0^ along with 


-Im 

TT 


Lio 1 -h 


m 


W 


ln(M' /m) 


W 


+ 4 Itifl + kj «(«'), (25) 

u' \ m J 


we hnd that the jet function to 0{as) is given by 




Un+p] 


U 


{mp + u' 


--In 1 + 


u 

rrir. 


d(u') + ( 1 + ) S(u') 


-I +4 

w 


\n{u' /rup) 


u' 


(26) 

The hrst line of (ESI) reduces to the result for the massless case in the limit rric — >■ 0 DI2I, 
while the second and third lines are unique to decay into charm quarks, and vanish for 
rric —>■ 0. 

We have also compared our calculation with the one-loop OPE result for b —> ciu in 
ES|.‘ For this purpose we re-expand the factorized expression for the hadronic tensor 
EH} in as- Using the notation of [2], the component kUi of the hadronic tensor, for 
instance, can be written as 

lUi 




HnS(u') + („') + 0[al ), 


(27) 


where the soft and jet contribution are given in ()2d|l and (l26|l . and the hard contribution 
reads j2] 


Hn = l + 


ttsCp 

An 


-AH -|- lOL — Alny 


2 In v , X i x 

^-4Li2 l-p ----12 28 

I — y 0 


^For the comparison with earlier calculations in 123 EH| see the detailed discussion in m 


























with y = nj^p/rub and L = In [ymb/y]. This has to be compared with the corresponding 
expression in [ 2 Bj 

= 7r ml Wu (29) 

where the resnlt for has to be expanded, nsing the SCET power connting ml/ml ~ 

-C 1. After some tedious, but straight-forward manipulations, we indeed 
hnd agreement with The comparison for the remaining components of the 

hadronic tensor is much easier because at order as they receive contributions from the 
hard functions only. Therefore the limit me/mb —0 in the corresponding expressions in 
pH] can be performed directly to recover the results for IT 4 5 from [2]. 


4 The partially integrated U spectrum 


From the results of the previous section we can derive any differential decay distribu¬ 
tion. We focus here on the U spectrum, because of its relation to the P+ spectrum from 
charmless decays. We pointed out in the last section that the hard and shape functions 
are unaffected by the presence of the charm-quark mass, and thus obey the same renor¬ 
malization group equations as in the charmless case. In the remainder of the paper, 
we will work with the renormalization-group improved formulas derived in j2]. After 
integrating over the lepton energy and neglecting higher-order terms in A, the doubly 
differential spectrum in the variables u and y is given by 

—^ / dujy‘^~''{Q-Ay)n{y)J{u-uj,mby,yi)S{uj,yi), (30) 

Tc dudy 


where —A < u < ymb —ml/ymb and me/mb < ?/ < 1. Note that after resummation all of 
the functions are to be evaluated at the intermediate scale pi ~ me- We have introduced 
the renormalization-group factors 

16 as{yi) 
a = — In —-—-, 

25 as[mb) 



and Vnimb, Pi), which resum logarithms between the hard and the jet scale. The exact 
form of Vh can be found in j2] , and the hard function Td can be derived from the functions 
Hij in the same reference. The total b ^ c rate to order as{mb) in the OPE is given by 


rc= IW 


cb 


(9f<\ \f(< 

\ 1927r3y 


+ 


CpCisimb) 

An 


25 

-- 2 .- 


me 

9 1 -2 


m 


b / J 


(32) 


At leading order in A we can set the phase-space factors /, g to unity, although higher- 
order corrections may be important numerically, as we shall discuss in Section ITHH 
It is useful to change from partonic to hadronic variables and dehne^ 


U = u + A, ci> = a;-l-A, S{uj) = S{uj) . 


(33) 


^Notice that uj is defined with the opposite sign in [2], whereas our convention for 5'(D) coincides. 
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(n_P) [GeV] 



(n+P) [GeV] 


Figure 2: Illustration of the shape-function region singled out by the cut on the variable 
U. The light-grey region shows the physical phase space M^/n+P < ri-P < n+P < 
Mb- The dark-grey part shows the shape-function region with A = 0.65 GeV and 
rric = 1.36 GeV. 


The relation between the hadronic momenta and the partonic momenta is given 
by n±P = n±p + A, where A = Mb — rub- This leads to 
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Tfl 

U = n_P -^ + C>(A^). (34) 

n+P 

To stay in the shape-function region, we need to restrict the phase-space integration to 
values of P ~ X^rrib- In order to preserve a close correspondence with the treatment 
of the B —>■ spectrum in [2], we introduce a cut U < A, with A being around 

600 MeV. The effect of this cut on the physical phase space in the variables P_ = n+P and 
P+ = n^P is illustrated in Fig. |21for typical values A = 0.65 GeV and rric = 1.36 GeV. 
The fraction of events with P < A is then given by 


Fc{A) = 


Fe(P < A) 

F. 


= / dA dy dU y‘^ - Ay)H{y) J{U - A, ynib) S{uj). (35) 


'nc Jo 

rtih 


A short calculation shows that the lower limit of the integration over y can be set to 
zero, up to terms of order {mc/mb)^~°'- After making this simplihcation, the integration 
limits are identical to those in 6 —>■ u decays. In fact, the integrals over the a* corrections 
from the hard function Pi and the hrst line of (ESI) are identical to the charmless case. 
We will give explicit results below. 

The new terms relevant to decay into charm quarks are contained in the last two 
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lines of fj2(i|l . After integration over U the result for these terms is 


/ dCjS{Cj) / dyy^-^^^Q-Ay) 


CpOisiyi) f 27r 


An 


- ln(j/A<2,) + 2\v?{y^co) 


1 + 2/A, 


+ ln(l + yAo,) + 4Li2(-2/Ai) L 


(36) 


where A;:^ = {A — u)mb/ml. The integrals over y can be evaluated in terms of the master 
integrals 


Gi{n,x) = dy 


y"- _ 2 Fi{l,n+ l;n + 2;-x) 

1 + xy n + 1 


G 2 (n,x) = / dy y"'\n(l + xy) = -( ln(l + x)— 

Jo n + 1 V 

^ , N Liof—x) G2in,x) 

G' 3 (n,x) = / dyy U 2 {-xy) = -+ 


+ 


Gi{n, x) 


n + 1 / n + 1 


n + 1 n + 1 ’ 
where the hypergeometric function 2 F 1 has a series expansion 

^ {di)k {0,2)k , N _T{ai + k) 

2Fi[ai,a2',a3;z) — 2 _^—Tyy- \o,i)k — 


(37) 

(38) 

(39) 


fc =0 


( 03 )?= kV 


r(ai 


(40) 


To express the hnal results in a compact way, we introduce 

6Gn(2 — a,X) — 4G„(3 — a, X) 


gn{a,X) = 

and make use of the functions dehned in [ 2 ] 

30 — 12a + 


/2(a) = 

T{a) = 


(6 — a) (4 — a) (3 — a) ’ 
2(6-a) 


T(a) 


/ 3 (a) = 


(41) 


2(138 - 90a TlSa^-a^) 
(6 — a)(4 — a)^(3 — a)^ 


(4 — a) (3 — a) 


(42) 


The hnal result can be written as Fc = Fu + Fm- The fraction is the result for b ^ u 
decays 


F„(A) = T(a)e 




duj >S(A, fjii) fu 


mfe(A — A) 


(43) 


f / \ 1 , Cpasimb) 

fu{x) = IH--- H{a) 


An 

GpOisiki) 

An 


[21n^x+ ( 4 / 2 (a) - 3 ) lnx+ (T-vr^ - 3/2(a) + 2 / 3 (a))] 
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An expression for H{a) can be found in j2]. The fraction is an additional piece unique 
to decay into charm quarks, which vanishes when rric —> 0. In the pole scheme, it is given 
by 


T™(A) 


T{a) 




duj >S(cD, p.j) fn 


mfc(A — ch) 


mi 


(44) 


fm{x) = 


CpOlg (/^i) 
47r 


2 In^ X + (4/2(0) — 1) In X 


27r^ 

+-^-/2(a) + 2/3(0) + 5(1(0, x) + 5(2(0, x) + 45(3(0, x) 


From the partially integrated spectrum -Fc(A) we can obtain the corresponding U spec¬ 
trum by differentiation, which results in 

F)1 

(45) 

where we have used integration by parts and A(0) = 0. 




fu 


mb{U — (jj) 


+ fi 


[mb{U 

m I 


m 


4.1 Change of charm-mass definition 

So far, our analysis has been performed with me dehned in the pole scheme. The effect 
of changing the charm-mass dehnition according to 


me —»■ me = me — 6m , 


where 6m ~ medslfii), is two-fold. First, the input value for the charm-quark mass 
in FA (A) is changed. Since the explicit charm-mass dependence in Fm is already an 
0{as) correction, this effect is formally of order Second, the jet function receives a 
perturbative correction proportional to 6m. It can be obtained from the tree-level jet 
function by taking into account the appropriate shift in the spectral variable, 

r/ N r /- N 2me6m ^ 

6(u — oj) ^ 6 (u — u) - 6 (u — uj), (46) 


where u = ri-p — ml/ripp is dehned using the new mass dehnition. Inserting the extra 
term into dSiD, one obtains an additional contribution to Fc(A), 


Fe{A) ^ Fe{A) - du dy dU - 4y) 


= Fe{A) - T(a + 1) S{A) . 

mb 


2me6m 

yrrib 


6'{U-u})S{u;) 


(47) 
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In order to see the scheme-independence of physical observables to a hxed order in 
one has to keep in mind that the relation between the hadronic momenta and the spectral 
variable U is also changed. Therefore, the result for Fc{A) in two different mass schemes 
should be compared at two different values of the cut-off parameter A, 


U = n_P 


mi r . 2fhr5m 


n+P 


< A ~ A -I- 


n+P 


such that Fc{A) in the new scheme reads 


Fe(A) = F„(A) + FrniA) - 6^^ T{a + 1) ^(A) . 


mb 


(48) 


(49) 


Expanding the upper limit A around A in the leading-order term in Fu{A) and neglecting 
terms of order we explicitly hnd the scheme-independence of our result, 

F,(A) = F,(A) + 0(a,2) . (50) 


Still, the convergence of the perturbative series at a given value of A might be rather 
different for different mass dehnitions. In addition to the pole scheme, we will consider 
two further examples, namely 

• the potential-subtracted (PS) scheme, where [20] 

’^c®(h/) =mc- /i/ + 0{al) (51) 

71 


with fif — 1 GeV, 

• the MS scheme, where 


mdfii) = me 


1 + 


CpO-s (/^i) 
dvr 


171 ^ 

3\n^-4] +0{ai) 

hi 


(52) 


4.2 Numerical predictions 

In this section we study the numerical predictions for Fc{A), taking into account mass- 
scheme and shape-function dependence. We start by summarizing the parameter values 
used in the subsequent analysis. The hard scale is hxed to the 5-quark mass, mb = 
4.65 MeV. The default value for the intermediate (jet) scale is /Xj = 1.5 GeV. We use 
the PS scheme as our default mass scheme, taking m^®(/i/ = 1 GeV) = 1.36 GeV. 
The charm-quark pole mass is taken as 1.65 GeV, and the MS mass at the jet scale as 
mdl^i) = 1-20 GeV. We use 2-loop running for with = 345 MeV, corresponding 

to asirrib) = 0.22 and as(/ij) = 0.37. 
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Figure 3: Predictions for partially integrated spectra in inclusive semi-leptonic b 
c decays: Left: NLO prediction for Fc{A) using the the default scenario S5 |2j in the PS 
scheme (solid line) compared to the MS scheme (long-dashed line) and the pole scheme 
(short-dashed line). Also plotted is the LO result (thick grey line). Right: NLO prediction 
for Fc{A) using the default scenario S5 in the PS scheme (solid line) compared to scenarios 
SI, S3, S7, S9 (dashed lines). 


For the numerical estimate we have to specify a model for the shape function, which 
we take from j2]^: 




^ r _ CpasiFi) f 
A dvr y 6 




(53) 


We use A = 0.685 GeV and b = 2.93 as our default (scenario “S5” in |2]). In Fig. El we 
compare the results for different mass schemes and different input shape functions as a 
function of the cut-off A. The following observations can be made: 


• For values of A ~ 600 MeV, the NLO corrections are large and positive. 

• Above some critical value Amax, the NLO corrections become so large that the 
fraction Fc exceeds 1, and therefore our result should not be trusted anymore. 

• The critical value A^ax amounts to about 480 MeV in the pole scheme, 700 MeV 
in the PS scheme, and 860 MeV in the MS scheme. 

• The model dependence from the input shape function amounts to an uncertainty 
of about 25%. 


4.3 Power corrections 

Our NLO calculation has been restricted to leading power in the 1 /mb expansion. Power 
corrections arise from two sources. First, one encounters new non-perturbative structure 

^ Notice that we have chopped off the radiative tail in S{uj), which does not contribute for the value 
of A that we are considering. 


14 















1.5 


power corrections 


1.4 
1.3 
1.2 
1.1 


0.2 0.4 0.6 0.8 


A [GeV] 


Figure 4: Effect of kinematic power corrections proportional to S{lo): The curve shows 
the NLO prediction for Fc{A) including the power corrections to the tree-level result in 
JSl, normalized to the leading-power result (j4.‘ll44|l . 


in the form of sub-leading shape functions. Second, there are kinematic power corrections 
proportional to ml/ml ~ and u ~ A^m^. The phase-space integration leads to 
logarithms ln(m^/m^) which can numerically enhance some of the power-suppressed 
terms.^ Whereas the estimate of sub-leading shape function effects is model dependent, 
the kinematic corrections multiplying the leading-order shape function can be calculated 
explicitly. We shall do this at tree level only, where we can use the results of to hnd 


r, du 



U-A 

mb 





+ 0{u\ A^) 


S{U) 


+ sub-leading shape functions . 


(54) 


The omitted terms are negligible in the portion of phase-space we are interested in. The 
numerical effect of the power corrections in (IH^ is plotted in Fig.|3J We see that Fc{A) is 
enhanced by about 20% at A = 0.65 GeV. Since we cannot control the remaining power 
corrections from sub-leading shape functions in a model-independent way, we consider 
this number as a rough estimate for the magnitude of the systematic uncertainties asso¬ 
ciated with power corrections. 


5 Relating b ^ c and b ^ u decays 

For the extraction of the CKM parameter |14fe| one would like to have a shape-function- 
independent relation between the B —>■ Xuii^t and B —>■ Xcii^e decay spectra. In what 
follows we focus on a relationship between the P+ spectrum in 6 —>■ n decays and the U 
spectrum in 6 ^ c decays. This relation can be obtained in a similar way as discussed 
for the comparison of P —> and B Xuii^e in ini. In the present case, this involves 

^These phase-space logarithms can be resummed using renormalization-group techniques, see m 
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Figure 5: Predictions for F„(A) using the weight function (1561) and the theoretical 6 —> c 
spectrum (ESI) on the basis of the shape-function model S5 in (p|). Solid line: PS scheme. 
Long-dashed line: MS scheme. Short-dashed line: Pole scheme. For comparison, we also 
show the direct computation of F„(A) from (l4(-il) (thick grey line). 


constructing a weight function W such that 


dP. 


dFu 

dP: 


F, 

F. 


dT 

dUW{A,U) 


__ 

dU "" ll/c 


dT 

dUW{A,U)^, (55) 


where we have used that F„/Fc = iKife/KbP to leading power in A. By measuring the 
partial decay rate F^i(P+ < A) mb —>■ u decays, as well as the dVc/dU spectrum in 6 —»• c 
decays, we can determine iKifel- The theoretical input is the weight function LF, which 
we can calculate from the results in ()44I49|1 : 


W{A,U) 


f mb{A-U) \ ^ CFas{f^i)T{a + l) ^ 

y 47r T(a) m, ^ ^ 

-|- 0{al) -|- power corrections, (56) 


in the PS scheme. We do not attempt to include power corrections here, but must be 
aware that they add a systematic uncertainty of at least 20%. 

At the moment, we do not have explicit experimental information on the U spectrum 
in P —> XciT^, so to illustrate how our method works we will have to rely on some 
theoretical input. In the following subsections we will consider two approaches: In 
the hrst we will use the theoretical prediction (ESD for the 6 —»• c spectrum to obtain the 
6 —>■ M spectrum from the weight-function analysis with dSEI)- In the second approach, we 
will construct a simple toy spectrum which takes into account possible charm-resonance 
effects. 


5.1 Numerical analysis using theoretical b ^ c spectrum 

In this sub-section we carry out the weight-function analysis using the theoretical b 
c spectrum from E3) as input. This will help us estimate some of the perturbative 
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Figure 6: Predictions for F„(A) using the weight function (1561) and the theoretical b ^ c 
spectrum (ESI) on the basis of the shape-function model S5 in (|53|) . The left plot shows 
the scale dependence (1 GeV< fii < 2.25 GeV) of the result using the PS scheme. The 
right plot shows the effect of varying by ±0.15 GeV around its default value. 


uncertainties inherent to our approach. We start by constructing the partially integrated 
6 —>■ M spectrum F„(A) from the theoretical 6 —>■ c spectrum (I45|l and the weight function 
W{A,U), using the shape-function model from In Fig. |S1 we compare the so- 

obtained results for F„(A) in the PS, pole, and MS schemes. We also show the result 
of the direct computation from (I4d|l . The difference between the curves is formally an 
0(ag) effect, and thus gives a rough measure of higher-order perturbative effects. We 
observe that in the pole scheme this effect is quite large for values of A below 700 MeV 
or so. At our reference point A = 650 MeV we obtain 

F„(0.65 GeV) = 0.71 from (HHl) and (IKH|l . PS scheme , 

Fu(0.65 GeV) = 0.62 from (HKll and (ESI), pole scheme , 

Fu(0.65 GeV) = 0.76 from (liH|) and (P^ . MS scheme , (57) 


compared to 

F„(0.65 GeV) = 0.79 from (EH|l , (58) 

from which we deduce a residual scheme dependence for F„(A) of about 10-15%. 

In Fig. ini we investigate the explicit /ij and rric dependence induced by the weight 
function (to isolate these effects, we £x rric and fii in the theoretical expression (j45j) for 
the b —>■ c spectrum). The charm-mass dependence is a small effect, less than 10% for 
reasonably large values of A. The dependence on the factorization scale fii is still sizeable, 
about 10-15%. The perturbative uncertainties related to the scheme and factorization- 
scale dependence could be resolved by calculating the corrections to the jet function 
in the massive case. 


5.2 Numerical analysis using a toy spectrum 

The purpose of this sub-section is to point out some aspects of the weight-function 
analysis that would be important when dealing with the physical 6 —> c spectrum. A 
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distinctive feature of this spectrum is that that the lowest-lying spin-symmetry doublet 
of charmed states D and D* already makes up about 80% of the semi-leptonic rate. For 
an on-shell D{D*) meson we formally have 


U, 


D{D*) 


^D{D* 


m: 


n. 


_P 


•C A {n+P ~ TUh) 


(59) 


Therefore, about 80% of the U spectrum is centered around “small” values of U (we put 
“small” in quotation marks, because numerically (M|, — ml)/mb ~ 350 MeV in the PS 
scheme). 

We will perform the weight-function analysis on a toy spectrum which takes this 
resonance structure into account. We construct this spectrum by assuming that the 
doubly-differential decay spectrum is concentrated along the D/D*-po\e and modulated 
by some function f{y), 


Tc d{n_P)dy \ n_P mt, j 


(60) 


where dyf{y) = 1, and Md = 1.975 GeV is a weighted average of the D and D* 
masses. We can derive the U spectrum (in a given mass scheme) from this model by 
taking 


n_P = U 




yrrib 


and performing the integral over y, which yields 


1 dTr 


M 


D 


m: 


r, du 


rubU^ 


f 


M 


D 


m: 


rribU 


9 U 


M 


D 


m" 


mb 


(61) 


For the following discussion, we use a simple parameterization 

^ T{2 + a + (5) _ p 

F(l + a)F(l + /3) ^ ^ 


(62) 


and £x the parameters a = 3.66 and (3 = —0.51 by requiring that Fc{3S.) at A = 650 MeV, 
and dVc/dU dXU = 550 MeV coincide with the theoretical expressions in the PS scheme.^ 
In Fig. [7| we compare the U spectrum from our toy model with the theoretical pre¬ 
diction in the PS scheme, using the shape-function model (j53|l as input. Fig. |H1 shows 
predictions for F„(A) obtained by applying the weight-function analysis to our toy spec¬ 
trum, as well as the theoretical curve obtained from (in. We see that at smaller values 
of A the sensitivity to the resonance structure and dependence on the mass scheme is 
sizeable. On the other hand, for larger values of A the resonance structure is washed 

®The reference point for J7(A) is sufficiently below the critical value A^ax = 700 MeV, and that for 
dVc/dU is sufficiently above the exclusive threshold C/min — 450 MeV. 
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Figure 7: The U spectrum in B —>■ Xcii'i decays: Theoretical prediction using the default 
model S5 for the shape function (dashed line) vs. phenomenological model (solid line) 
assuming the dominance of a single D-meson pole. The model parameters are adjusted 
to reproduce the value of the spectrum (left) at U = 550 MeV as well as the integrated 
spectrum (right) at A = 650 MeV. 
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Figure 8: Predictions for the partial rate F„(A) in .B —> from a toy spectrum in 

B XciH'i using the NLO weight function. Comparison of PS scheme (solid line), pole 
scheme (dotted), and MS scheme (dashed) with the theoretical result (jdhl) using scenario 
S5 (thick grey line). 


out, and the predictions obtained in different mass schemes converge. The sensitivity 
to the resonance structure at moderate values of A means that the phenomenologically 
acceptable window for A in the shape-function approach is smaller than in the b ^ u 
case, where the contributions from the charmless ground states vr, rj, p and uj add up 
to only about 25% of the total semi-leptonic b ^ u rate and are moreover centered at 
values of P+ not much larger than 100 MeV. 

These observations have important implications for extracting \Vub\ by relating par¬ 
tially integrated b ^ u and b ^ c decay spectra. To apply our results to 6 —> c decays, 
it is crucial that the cut-off parameter A be sufficiently large to avoid sensitivity to the 
shape of the spectrum in the resonance region. To apply them to b ^ u decays, A must 
be small enough to suppress the charm background, which sets in at A ~ 650 MeV. 
Balancing between the two cases restricts the cut-off parameter A to a rather small 
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window. 

We also observe that the weight-function analysis with our toy model systematically 
underestimates the result for F„(A) compared to the “true” result For instance, 

at our reference point A = 0.65 GeV, we have 

F„(0.65 GeV) = 0.79 from , (63) 

F„(0.65 GeV) = 0.55 from toy model and PS scheme . (64) 

This is due at least in part to the crudeness of our model, which completely ignores the 
non-negligible continuum contribution. While we could refine our model to take this into 
account, we think that such hne-tuning is best resolved by experimental input. 

6 Conclusions 

We analyzed perturbative corrections to i? —>■ Xcivu decays using the power counting 

~ a/Aqcd^ 6 for the charm-quark mass. This treatment implies that a certain class 
of partially integrated b ^ c decay spectra is sensitive to the non-perturbative shape- 
function effects familiar from B —> Xuii^i decays. With the aid of soft-collinear effective 
theory, we showed that the one-loop corrections to such decay spectra can be written 
as a convolution of hard, jet, and shape functions. The hard and shape functions are 
identical to those found in the factorization formula for B —> X^ii^i decays, but the jet 
function depends explicitly on rric and hence receives non-trivial corrections unique to 
decay into charm quarks. We calculated these corrections at NLO in perturbation theory 
and at leading order in the l/rrib expansion, and derived a shape-function independent 
relation between partially integrated B —> Xcii^e and B —>■ Xuii^i decay spectra. This 
relation can be used to determine | 14 fe|- 

Numerical studies raised some issues related to this treatment. First, the portion 
of phase-space where the shape-function approach is valid is somewhat smaller in i? —> 
Xcii^i decays than in the charmless case. Second, although the results are formally 
independent of the renormalization scheme used to define the charm-quark mass, the 
numerical dependence on the mass scheme is signihcant. Finally, the structure of power 
corrections is slightly more complicated than in the charmless case, since one encounters 
not only sub-leading shape functions, but also kinematic power corrections. Some of the 
power corrections are enhanced by large logarithms ln(m^/mj). 

Our study may help improve the understanding of inclusive B decays in the shape- 
function region. On the one hand, it provides additional information for the extraction 
of \Vub\- On the other hand, it may offer an additional testing ground for theoretical 
methods based on factorization and soft-collinear effective theory. To explore these ideas 
further would require experimental information on the partially integrated B Xcii^e 
decay spectrum used in our analysis. 
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